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ABSTRACT 

Radiation pressure-dominated accretion disks are predicted to exhibit strong density inhomogeneities 
on scales much smaller than the disk scale height, due to the nonlinear development of photon bubble 
instability. Radiation would escape from such a "leaky" disk at a rate higher than that predicted by 
standard accretion disk theory. The disk scale height is then smaller than that of a similar disk without 
small-scale inhomogeneities, and the disk can remain geometrically thin even as the flux approaches and 
exceeds the Eddington limit. An idealized one-zone model for disks with radiation-driven inhomogeneities 
suggests that the escaping flux could exceed I/^ by a factor of up to 10 — 100, depending on the mass 
of the central object. Such luminous disks would develop strong mass loss, but the resulting decrease 
in accretion rate would not necessarily prevent the luminosity from exceeding Le- We suggest that 
the observed "ultraluminous X-ray sources" are actually thin, super-Eddington accretion disks orbiting 
stellar-mass black holes, and need not indicate the existence of a class of intermediate-mass black holes. 

Subject headings: accretion: accretion, accretion disks - black hole physics - hydrodynamics - MHD - 
X-rays: binaries 



1. INTRODUCTION 

Magnetized, radiation pressure-dominated atmospheres 
tend to become highly inhomogeneous on scales much 
smaller than the radiation pressure scale height, due to 
the nonlinear development of "photon bubble instability" 
(Arons 1992; Gammie 1998). In this mechanism, grav- 
ity squeezes the radiation out of slightly overdense re- 
gions, while magnetic tension enhances the overdensity by 
preventing the gas from spreading sideways. Simultane- 
ously, slightly underdense regions grow more tenuous as 
radiation force accelerates them; hence the density con- 
trast increases. In a nonlinear analogue of this instability 
(Begelman 2001; hereafter. Paper I), gas cycles between 
narrow regions of very high density and much broader re- 
gions of very low density. Gravity and gas pressure domi- 
nate the force balance in the dense regions, which resem- 
ble gas pressure-supported exponential atmospheres, while 
the low density regions resemble radiation pressure-driven 
winds constrained to slide along magnetic field lines. To 
complete the cycle, low density matter rejoins the dense 
phase by passing through an isothermal, gas pressure- 
dominated shock, while gas moves from the high-density 
phase to the tenuous phase more gradually, via a slow mag- 
netosonic critical point. Thus, the inhomogeneities repre- 
sent nonlinear wave patterns rather than discrete, long- 
lived filaments. 

Thanks to radiative diffusion, the presence of large den- 
sity contrasts has little effect on the radiation pressure gra- 
dient. Under such conditions the Eddington limit can be 
circumvented. In the diffusion limit, the radiation flux is 
inversely proportional to the local density. Radiation flows 
readily through tenuous regions while avoiding regions of 
high density. If the low- and high-density regions are cou- 
pled dynamically, the upward radiation force in the ten- 
uous zones can balance the downward gravitational force 



acting on the dense regions, and global dynamical equilib- 
rium is maintained. In the steady-state, nonlinear plane- 
wave solutions presented in Paper I, a combination of ram 
pressure and magnetic tension forces maintains this cou- 
pling. If most of the mass resides in the dense phase while 
most of the volume is filled with tenuous gas, the total flux 
through the system can exceed the Eddington limit by a 
large factor (Shaviv 1998, 2000). 

This effect should be relevant in luminous accretion 
disks, the inner regions of which are radiation pressure- 
dominated (Shakura & Sunyaev 1973). If there are no 
small-scale inhomogeneities, as assumed in the standard 
theory, the vertical scale height /i at a given radius is pro- 
portional to the local radiation flux. If the luminosity ap- 
proaches the Eddington limit for the central mass, such a 
disk must inflate to a thickness of order its radius, proba- 
bly triggering substantial mass loss. But in the presence of 
strong density inhomogeneities the disk becomes "leaky," 
and more radiation can escape without inflating the disk. 
One can even imagine a scenario in which the luminos- 
ity exceeds the Eddington limit, while the disk remains 
geometrically thin. 

In this Letter we analyze the possible effects of nonlinear 
photon bubble instability on thin accretion disks. In § 2 
we summarize some generic features of the resulting inho- 
mogeneities in radiation-dominated atmospheres. In § 3 
wc consider their effects on the structure and luminosity 
of accretion disks, and in § 4 we discuss the results and 
propose possible applications to astrophysical systems. In 
particular, we suggest that the observed "ultraluminous 
X-ray sources" are accretion disks radiating at above the 
Eddington limit for their central masses. 
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Consider a plane-parallel radiation-dominated atmo- 
sphere with column density E and scale height h in a, grav- 
itational field g{h). The radiation flux passing through 
the atmosphere is F and the Eddington flux is F^ih) ~ 
cg{h)/ K, where k is the opacity. The mean mass density is 
p = E/2/i. If most of the volume of the atmosphere is filled 
with tenuous gas with density <C p, then the condition 
for overall dynamical equilibrium is (Shaviv 1998) 
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where Pr is the radiation pressure. Since gas cycles be- 
tween low and high densities via isothermal shocks, there 
is an additional, local pressure balance condition between 
the ram pressure of the tenuous phase, p-v'^ , and the post- 
shock pressure, p+c^, in the dense phase, where Cg is the 
thermal sound speed in the gas. Downstream of the shock, 
matter accelerates along each magnetic flux tube (which 
is assumed to be approximately straight on scales of or- 
der the gas pressure scale height, c^/g <C h), and passes 
through a sonic point where Vs = Cg (Paper I) . This implies 
a constant mass flux rhs — PsCg, where ps is the density 
at the sonic point. If radiative transfer obeys the diffusion 
law with F (x p^^, then ps — p. We therefore have 
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Thus, p is the geometric mean between p_ and p+. 

We solve the two pressure balance equations for 
in terms of F and the other parameters. These results 
are verified by detailed calculations of steady-state, plane- 
parallel inhomogeneities (Paper I). This yields the "Ed- 
dington enhancement factor," £ = F/Fe ^ P+/p, giving 
the factor by which the radiation flux through the atmo- 
sphere exceeds that through a homogeneous atmosphere 
with the same column density and scale height. 

The above analysis applies only if the standard radia- 
tive diffusion equation is valid locally throughout the flow, 
i.e., if the gas is optically thick across any density scale 
length along the flow, p^K/|Vp| > 1. Where this condition 
is violated the radiation field does not "see" the change in 
density and the flux does not adjust to local conditions. 
The local diffusion condition is first violated near the min- 
imum density, p_ (Paper I), placing a lower limit on the 
critical mass flux, ^ g/cgK. In this limit, the critical 
point occurs at a density larger than p. The Eddington 
enhancement factor is therefore 
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To complete the model we specify p+ in terms of other 
atmospheric parameters. While the nonlinear theory in its 
current state does not predict the level of inhomogeneity, 
it clearly indicates an upper limit. If the maximum gas 
pressure p+Cg exceeds the magnetic pressure, the fleld will 
buckle and the crucial dynamical coupling will be lost. We 
conjecture that this limit also represents the natural satu- 
ration level of inhomogeneous structure. We therefore set 
p+Cg ~ £,Pm, where Pm is the magnetic pressure and ^ < 1 
is a parameter reflecting the "inefficiency" of the process 



by which inhomogeneities grow (e.g., due to the geometry 
and/or time-dependence of the magnetic field, etc.). The 
value of ^ is not known, but we guess that it might lie in 
the range ~ 0.01 — 0.1. In anticipation of applying our 
model to accretion disks, we set pm ~ apr where a might 
also lie in the range ^ 0.01 — 0.1. Therefore, from eq. (||) 
we have 
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where pg = pCg <C Pr is the mean gas pressure. The cor- 
responding spacing between dense regions — the "wave- 
length" of the inhomogeneities — is given by A 
^ahmin(l,pgK/g), which is comfortably smaller than the 
radiation pressure scale height. Note that Cg can be esti- 
mated by assuming that the gas is in LTE with the radi- 
ation field. The correct a to use here corresponds to the 
total magnetic pressure, rather than the BrB^ stress asso- 
ciated with the transport of angular momentum. Simula- 
tions suggest that the latter might be ^ 10 times smaller 
than Pm, which is dominated by B^/Stt (Brandenburg et 
al. 1995; Stone et al. 1996). In applying this model to ac- 
cretion disks, we correct for the discrepancy by choosing 
a larger fiducial value of ^, ^ ~ 0.1. 

3. APPLICATION TO ACCRETION DISKS 

Now consider the atmosphere in the context of a thin, 
Keplerian accretion disk at radius r with accretion rate 
M. We will primarily consider a one-zone model in which 
the disk is characterized by a single mean density and 
scale height; later we will consider the effects of density 
stratification. We normalize r to the gravitational radius 
of a central mass M, M to the Eddington accretion rate 
Me — 4:TtGM/ Kc, and h to the radius: 



M , h 
m = —. — ; d = —. 
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The local flux due to accretion is then 
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F=- ^ — V = - x-^mD, (6) 
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where P = 1 — (xi-a^/xY/'^ for the standard assumption of 
zero torque at the innermost stable orbit, Xin (Shakura & 
Sunyaev 1973). The Eddington enhancement factor is 
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In a standard radiation-dominated accretion disk 5 adjusts 
so that t — 1, but this will not be the case in the presence 
of small-scale inhomogeneities. For an a— model viscosity, 
the dissipated flux is given by 

F « -aY.n^'e = l^ai:5^x-^'\ (8) 
2 2GM ^ ' 

Equating the accretion flux and the dissipated flux we ob- 
tain 
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The radiation and gas pressure are calculated assum- 
ing vertical hydrostatic equilibrium and LTE, respectively. 
Setting K = 0.4 cm^ g""'^ for electron scattering in ionized 
hydrogen, and letting m = M/Mq, we have 
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Combining these formulae with equations (Q) and (Q), 
obtain 



we 
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where we have normalized ^ to 0.1 and a to 0.01. A similar 
equation for i can be derived from eq. (p^), using eq. (^. 
Inhomogeneities affect the disk structure only if ^ > 1, 
implying 
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where we have set I? 1 for Xinhom S> x-^. 

Now consider the maximum possible mass flux that can 
be carried by the disk. Since inhomogeneities are ineffec- 
tive outside Xinhom, we have mmax(iCinhom) ~ 2a:;inhom/3, 

corresponding to a thick disk [^(xinhom) 1] with a lo- 
cal radiation flux equal to the Eddington limit. Any 
excess matter supplied at larger x presumably will be 
driven away by radiation pressure before reaching Xinhom- 
For m close to rhmax and ^ < 0.1, local radiative dif- 
fusion should apply through most of the inhomogeneous 
region, and certainly in the inner regions where most 
energy is released. We therefore have mmax(a;inhom) 

580^^'^]^^a?.2?7i^/^. By setting S{x) ~ 1 at x < xinhom, 
we find that the carrying capacity of the disk decreases 
with decreasing r, rhmax oc x^^'^/'D, down to the ra- 
dius X 4xin at which x^^'^/D has a local minimum. 
Thus, the mass flux reaching the central object may be 
much smaller than that entering the inhomogeneous re- 
gion, min,max 4(a;in/a;inhom)i^^n2max(a;inhom)- Neverthe- 
less, the luminosity can exceed Le- 

According to the one-zone model, the maximum possible 
disk luminosity satisfies 
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where e is the accretion efficiency and the inequality re- 
flects the fact that m increases with r. However, mass 
loss from the disk's upper layers may limit the maximum 
attainable luminosity. Such a wind must arise, in a super- 
Eddington disk, from strata where the mean optical depth 
is so small that density inhomogeneities provide insuffi- 
cient modulation of the flux. If the radiation becomes 



"trapped" in the accelerating outflow it will drive mass loss 
at a rate comparable to or exceeding the inflow rate, pos- 
sibly regulating the accretion luminosity to a value ~ L^. 
But if the gas remains sufficiently inhomogeneous in the 
acceleration zone, radiative losses due to diffusion could 
reduce the energy supply available for unbinding the gas. 
The mass loss rate would then be lower (Shaviv 2001). 
Mass loss could be further inhibited if dynamo processes 
in the disk maintain a sufficiently dense array of closed 
magnetic flux loops in a magnetically dominated corona 
(Miller & Stone 2000). If £ < Pmn/g, then gas propelled 
upward into closed field lines anchored in the disk can be- 
come trapped by the magnetic tension. Once its optical 
depth reaches ~ ^, a trapped blob begins to be weighed 
down by gravity, and if r reaches ^ PmK/ g before the field 
opens up, the gas will overcome the tension and drop back 
onto the disk. 



4. DISCUSSION 

We have argued that the radiation-dominated, inner re- 
gions of accretion disks may be able to approach or even 
exceed the Eddington limit, without becoming geometri- 
cally thick or suffering catastrophic mass loss. The disks 
become highly inhomogeneous on small scales, and corre- 
spondingly porous to radiation leakage, due to the non- 
linear development of photon bubble instability (Gammie 
1998; Paper I). The possible enhancement in radiation flux 
is limited by the strength of the magnetic field, which in 
turn is related to viscous dissipation and angular momen- 
tum transport via the magnetorotational instability (MRI, 
which operates concurrently with photon bubble instabil- 
ity [Blaes & Socrates 2001; Turner, Stone, & Sano 2001]). 
Greater porosity requires larger density constrasts, which 
place greater stresses on the magnetic field. Beyond a cer- 
tain flux, the magnetic tension will be overwhelmed and 
the field will buckle. Simple one-zone models suggest that 
accretion disks around stellar-mass black holes could radi- 
ate as much as ^ 10 times the Eddington limit. A weak 
mass dependence implies that disks around supcrmassive 
black holes could exceed the Eddington limit by an even 
larger factor, perhaps ~ 100 or more. 

Although the magnetic configurations in accretion disks 
are highly dynamical and do not resemble the quiescent, 
uniform fields assumed in our analytic models, this time- 
dependence need not impede the development of inhomo- 
geneities. The linear growth time for photon bubble in- 
stability is ~ (A/g)"'^/^, the free-fall time across the wave- 
length (Gammie 1998). This is shorter than the orbital 
time, which determines the evolution timescale of mag- 
netic structures due to MRI (Blaes & Socrates 2001). 

Catastrophic mass loss from the disk's upper layers 
poses the most serious threat to the production of super- 
Eddington luminosities. If the luminosity is effectively 
trapped in the outflowing gas it could deplete the accre- 
tion flow, regulating the maximum luminosity to Le- 
But if inhomogeneities persist into the acceleration zone 
the mass loss could be lower, leaving the accretion flow 
largely intact. A magnetically structured corona with a 
sufficient number of closed field lines could also inhibit 
mass loss by trapping and recycling some of the expelled 
gas. Both speculations will have to be checked by detailed 
calculations, including numerical simulations. 
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Fluxes ^ IOLe could explain the existence of "ultralu- 
minous X-ray sources" (ULXs: Fabbiano 1989; Colbert & 
Mushotzky 1999; Makishima et al. 2000; Fabbiano, Zezas, 
& Murray 2001). The high luminosities of these sources 
have led to suggestions that they represent a new class of 
intermediate-mass black holes. However, the observed lu- 
minosities are also explainable in terms super-Eddington 
fluxes from disks around black holes of a few solar masses. 
This is consistent with anomalies reported by Makishima 
et al. (2000), who performed "multicolor blackbody disk" 
fits to ULX spectra and found that the best-fitting models 
were considerably hotter and more compact than expected 
for ~ lOOAf0 black holes. King et al. (2001) have also pro- 
posed that ULXs are sources producing super-Eddington 
fluxes, but they posit that the intense radiation is beamed 
within a small solid angle. In contrast to our proposal, 
which predicts a roughly isotropic flux, they require a large 
population of misdirected sources. However, we agree with 
them that these systems should not be very long-lived, 
since an object accreting at ~ 10 times the Eddington 
limit would double its mass in a few million years. 

We make two other predictions based on the analysis 
presented in this paper. First, inhomogeneous accretion 
disks radiating above the Eddington limit should possess 
"cold" winds driven by continuum radiation pressure. In 
contrast to thermally driven outflows, these winds should 
be at or below the Compton temperature of the accre- 
tion disk radiation, and well below the virial tempera- 
ture. They should have substantial optical depths and may 
therefore modify the underlying disk radiation. Moreover, 
if "hot" flares (energized by magnetic reconnection, for 
example) percolate through such a wind, Compton down- 
scattering of the hard radiation could produce a prominent 
"soft excess" . It may be difficult to distinguish such winds 



from large-scale hydromagnetic flows (Blandford & Payne 
1982; Konigl & Kartje 1994), if the latter also possess large 
optical depths. 

Second, there may be a population of "ultraluminous" 
AGNs, i.e., supermassive black holes radiating at well 
above their Eddington limits. Candidates could be iden- 
tified by looking for discrepancies between the black hole 
mass as estimated from the properties of the host galaxy 
(Ferrarese & Merritt 2000; Gebhardt et al. 2000) or by re- 
verberation mapping and line widths, and the lower bound 
on the mass according to the Eddington limit. An epoch 
of super-Eddington accretion at high z could lead to the 
rapid growth of supermassive black holes by accretion. 

Finally, we note that SS 433 may be an example of an 
object which is accreting at a rate that far exceeds even 
the upper limit attainable in a "leaky" disk (King, Taam, 
& Begelman 2000). Most of the matter initially in the 
accretion flow is therefore blown away before it reaches 
the vicinity of the central object. Nevertheless, our anal- 
ysis suggests that the mass flux reaching the central ob- 
ject could exceed the Eddington limit by a substantial fac- 
tor. Much of the energy released by the accreting matter 
would be available for driving an outflow and perhaps goes 
into powering the 0.26c jets, which are believed to carry a 
super-Eddington energy flux (Begelman et al. 1980). We 
presumably do not see a highly super-Eddington radiative 
flux because of the large optical depth of intervening mat- 
ter. 
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